A Chandra ACIS Observation of the X-ray Luminous SN 1988Z 
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ABSTRACT 

SN 1988Z is the most luminous X-ray-cmitting supernova, initially detected in 1995 using 
the ROSATBRl with a luminosity of -8x10''" erg s^^ (Fabian & Terlevich 1996). Its high 
luminosity was ascribed to expansion of the blast wave into an especially dense circumstellar 
medium. In this paper, we describe a recent observation of SN 1988Z using the ACIS detector 
on Chandra. We readily detect SN 1988Z, obtaining '^30 net counts which corresponds to a 
0.2-2.0 keV luminosity of —3.2x10^^ erg s^^. The calculated quantiles for the extracted counts 
allow a broad range of temperatures, but require a temperature hotter than 5 keV if there is no 
intrinsic absorption. The long term light curve (1995-2005) declines as t~^-^='="-^. This is one of 
the steepest X-ray light curves. The X-ray luminosity indicates that the emitting region has a 
high density (>10^ cm~^) and that the density profile is not consistent with a constant mass loss 
stellar wind during the —5,000 years before the explosion. If the circumstellar medium is due 
to progenitor mass loss, then the mass loss rate is extremely high (— lO^'^ Mq yr^^(v,uj/10 km 
s~')). The X-ray results are compared with the predictions of models of SN1988Z. 



Subject headings: supernovae: individual (SN 1988Z) 



Introduction 



The detection of X-rays from a supernova is a 
rare event. While several mechanisms might pro- 
duce X-ray emission - a thermal flash associated 
with breakout of the shock from the stellar sur- 
face, Compton scattered 7-rays from metals syn- 
thesized in the explosion, shocked-heated circum- 
stellar matter - only the latter has been detected 
more than once. Down-scattered 7-rays have been 
detected only from SN 1987A between 6 and 12 
months after the explosion. Strong X-ray emis- 
sion from shock-heated circumstellar matter oc- 
curs only when a relatively large amount of mate- 
rial is located in close proximity to the progenitor. 
As this material is generally the product of the 
pre-supernova stellar mass loss, measuring the X- 
ray light curve reveals the mass loss history. As the 
blast wave propagates outward considerably faster 
than the wind, the supernova light curve over the 
course of a few years reveals the mass loss history 
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of the progenitor's final several thousand years. 

The total list of supernovae detected in X-rays 
consists of —24 objects, and is slowly growing 
as a result of Chandra, XMM-Newton and now 
Swift observations. ' Long term light curves exist 
for about half a dozen of these (Inimler & Kuntz 
2005). These light curves show a range of behav- 
ior, from nearly flat (e.g., SN 1979C and SN 1978K 
since its 1991 discovery), to declining at various 
rates, to more complex behavior (e.g., SN 1993 J - 
Zimmermann & Aschenbach (2003)). These light 
curves are decidedly incomplete, however, as data 
sets consisting of more than four observations exist 
for only four SNe (SN 1978K, SN 1993J, SN 1998S, 
SN 1999em) and only six supernovae have been fol- 
lowed beyond a few years. In this paper we add 
SN 1988Z to this select group. 

SN 1988Z is one of a small number of extremely 
luminous, X-ray emitting supernovae (Schlegel 
1995; Immler & Lewin 2003; Schlegel 2005) 
and is the most distant confirmed X-ray emit- 
ter. SN 1988Z was discovered on 1988 Dec 12 
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by two independent observers, C. PoUas (PoUas 
1988) and G. Candeo (Cappellaro, Turatto, & 
Candeo 1988), near the galaxy MCG +03-28-022 
(= Zw 095-049). It was classified as a Type II 
by Heathcote, Cowley, & Hartwick (1988) from 
an optical spectrum. Filippenko (1989) confirmed 
the classification, noting the spectral peculiari- 
ties of very blue color and lack of P Cyg profiles. 
Schlegel (1990) included it as a member of a new 
subclass, the Type Iln supernovae. Weiler, Pana- 
gia, & Sramek (1990) reported the first radio 
detection of SN 1988Z using the Very Large Ar- 
ray; subsequent detailed early radio observations 
were reported in Van Dyk et al. (1993) and re- 
cent radio observations in Williams et al. (2002). 
Aretxaga et al. (1999) amassed the observations 
of SN 1988Z available as of late 1998 and con- 
structed the first spectral energy distribution for 
any Type Iln supernova. 

SN 1988Z's importance arises from its extreme 
distance, ^93 Mpc (assuming Ho=72 km s^^ 
Mpc~^). ^. At a discovery luminosity of '-^8x10'*^ 
ergs s~^, it is also the most luminous. (This lumi- 
nosity was revised down to (3±1.6)xl0^° ergs s^^ 
by Aretxaga et al. (1999)). More importantly, 
SN 1988Z has been proposed by Fabian & Ter- 
levich (1996) as the prototype "compact super- 
nova," an event in a high density medium emit- 
ting sufRcient luminosity to power the broad-line 
regions of some active galactic nuclei (e. g., Ter- 
levich et al. (1992)). While this idea has fallen 
out of favor, the existence of SN 1988Z suggests 
that extremely luminous supernovae occur via the 
proposed mechanism, and are interesting objects 
of study in their own right. A second interpre- 
tation suggests that the high luminosity arises 
from shocks driven into dense clumps formed 
by an unstable stellar wind (Chugai & Danziger 
(1994);Winiams et al. (2002)). Either mechanism 
requires a substantial circumstellar envelope. 

In this paper, we present the result of a Chandra 
ACIS observation of SN 1988Z. C/iandra 's sensitiv- 
ity allows us to clearly detect the object a decade 
after its initial X-ray detection. While limited 
spectral information is available from the observa- 
tion, the inferred flux and the 10-year light curve 
lead to interesting inferences about SN 1988Z and 



^A marginal detection of SN 2002hi, at an estimated dis- 
tance of 180 Mpc, is reported by Pooley & Lewin (2003) 
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its environment. 

Throughout this paper we adopt a distance to 
SN 1988Z of 93 Mpc, calculated from the redshift 
of 6670 km s^^ (Barbon, Cappellaro, & Turatto 
1989) and using the value for Hq of 72 km s^^ 
Mpc~^ derived from both WMAP (Spergel et al. 
2003) and the Hubble Space Telescope (Freedman 
et al. 2001) . Note that this is different from 
the distance assumed in previous papers describ- 
ing X-ray and radio observations (Fabian & Ter- 
levich (1996); Aretxaga et al. (1999); Wilhams et 
al. (2002); Table 1). 



2. Observation 



Chandra observed SN 1988Z on 2004 June 29 
(MJD 53186) for -50 ksec with the ACIS detector 
(Garmire et al. 2003). The aimpoint of the de- 
tector fell on the back-illuminated S3 CCD. The 
observation was cut short by ~2 ksec by increasing 
electron contamination following a solar flare and 
during approach to the radiation belts. We ex- 
tracted a light curve of the background, masking 
out all bright sources detectable by eye, to search 
for any signatures of flaring events known to af- 
fect back-illuminated CCDs (Chandra Proposers' 
Guide 2002). The background rate increased dra- 
matically near the end of the observation, reflect- 
ing the increasing radiation background from the 
solar flare. Removal of the increased background 
interval yielded a net exposure time of 44.62 ksec. 



Results 



Figure 1 shows the optical field containing SN 1988Z 
and its host galaxy MCG +03-28-022. Circles in- 
dicate the X-ray sources detected in the 0.3-10.0 
keV band. A total of four sources were detected 
within 2 arc minutes of the nucleus. The brightest 
of these is located within 0".5 of the radio position 
of SN 1988Z (Wilhams et al. 2002). This source 
is marked in Figure 1 with the white circle. It is 
the brightest source among the four by a factor of 
two. Note that the nucleus of the host galaxy is 
not detected as an X-ray source. 

We extracted counts from SN 1988Z using an aper- 
ture of 4" centered on the source; this aperture 
encloses >98% of the point spread function of an 



on-axis point source (Chandra Proposers' Guide 
2002). We also extracted counts in several bands; 
for example, to match the bandpass of the ROSAT 
HRI (Aretxaga et al. 1999), we used the 0.2-2 keV 
band. There is a larger uncertainty (and fewer net 
counts) in the 0.2-2 keV band than in the 0.3-2 keV 
band because of the higher background below 0.3 
keV (Chandra Proposers' Guide 2002). 

The background was extracted from a circle of ra- 
dius 20" centered ~30" to the southeast of the 
source circle to avoid any possible emission from 
the host galaxy. The net total-band count rate 
is (8.4±1.7)xl0-'* counts s^^. The rates for all 
extracted bands are listed in Table 2. At these 
low count rates, the effects due to event pileup are 
negligible (Davis 2001). 



3.1. Spectral Constraints 

Quantiles have been demonstrated to be more ro- 
bust than hardness ratios for estimating spectral 
properties of sources with very limited numbers of 
counts because they do not generate upper limits 
(Hong, Schlegel, & Grindlay 2004). For the ACIS 
observation, we extracted the photon energy val- 
ues for the events in the source and background 
apertures and generated the 25% and 75% quan- 
tile values. 

Figure 2 shows the result of quantile analysis for 
a coronal ionization equilibrium plasma (MeKaL) 
grid. The MeKal grid has the depicted shape be- 
cause of the changing presence of line emission as 
the temperature (primarily) is varied. The quan- 
tile analysis yields a preferred temperature of ^2- 
3 keV, but the 90% uncertainty covers a very 
broad range, leaving the temperature high, but 
essentially undefined. The preferred column den- 
sity is 5x10^^ cm^^, with a range between 1x10^^ 
cm~^ and 2x10^^ cm~^. The broad allowed tem- 
perature range can be understood from the rel- 
ative instability of quantiles using the MeKaL 
model caused by the emergence and decline of var- 
ious lines as the temperature varies. If we exter- 
nally constrain the value of the column density by 
using the Galactic value of 1.4x10^° cm~^, then 
the temperature is constrained to be hotter than 
~5 keV. 



3.2. Flux and Light Curve 

Converting the measured count rates into physi- 
cal units requires estimates of distance and col- 
umn density plus an assumption about the emis- 
sion mechanism. We use different values for all 
three from those used in previous papers (Fabian 
& Terlevich (1996), Aretxaga et al. (1999)), as 
justified below. 

We estimated the column density toward SN 1988Z 
using the dust maps of Schlegel, Finkbeiner, & 
Davis (1998) and the Eb_v-Nh relation of Pre- 
dehl & Schmitt (1995). Together these yield a 
value of 1.4x10^° cm^^, lower than the 2.5x10^" 
cm~^ from the lower angular resolution maps of 
Burstein & Heiles (1982) used by Aretxaga et al. 
(1999). Aretxaga et al. (1999) point out that the 
X-ray column density in the immediate vicinity 
of the supernova could be much higher if the pre- 
ionized material enshrouding the supenova is only 
partially ionized, as the partially ionized material 
will absorb a substantial fraction of the soft X-ray 
flux. They estimate that the column density could 
easily be as high as 10^^ cm^^. The absence of 
absorption in the radio data suggest this is not 
the case (Williams et al. 2002). Nevertheless 
it should be noted that the soft X-ray luminosity 
would be substantially underestimated if we adopt 
the Galactic column density while the supernova 
instead has a high intrinsic column. 

A supernova shock encountering dense circumstel- 
lar material or ejecta generally produces a line-rich 
X-ray spectrum from atomic transitions within 
nearly stripped ions. Given the inferred density 
of the shocked medium around SN 1988Z, it is ex- 
pected to rapidly attain collisional ionization equi- 
librium (Chevaher & Fransson 1994, 2004). We 
therefore use the MeKaL model for estimating flux 
and luminosity (Liedahl et al. 1992; Mewe et al. 
1985). This is also different from previous publica- 
tions, which used thermal Bremsstrahlung (Fabian 
& Terlevich 1996; Aretxaga et al. 1999) (Ta- 
ble 1). If SN 1988Z emits copious hard X-rays, 
then it is possible that the entire surrounding 
medium has become fully ionized, in which case 
Bremsstrahlung emission would dominate. Given 
the absence of spectral information, either spectral 
model can be used equally well to infer relative 
fluxes among the available X-ray observations. 



The quantile analysis indicates that if we exter- 



nally constrain the value of the column density by 
using the Galactic value of 1.4x10^° cni~^, then 
the temperature is constrained to be hotter than 
~5 keV. To maintain consistency with the quantile 
results we assume for our flux estimates a 5 keV 
plasma with Galactic column density. Flux analy- 
sis of the X-ray emission from other remnants with 
higher flux has generally used a temperature of 0.8 
keV, inferred from spectral fitting (Schlegel 2005; 
Immler et al. 2005). This has been intepreted 
as the temperature of the reverse shock. It should 
be noted that the emissivity of a plasma in the 
0.2-2.0 keV band varies slowly over the temper- 
ature range from 0.5 to 10 keV. The emissivity 
for a 0.8 keV plasma is approximately 25 percent 
higher than a 5 keV plasma. The acknowledged 
inconsistency with other analyses resulting from 
the use here of the higher temperature is however 
comparable with or smaller than other uncertain- 
ties, such as possible temperature evolution be- 
tween the RO SAT and Chandra observations. 

The recession velocity of MCG +03-28-022 is 6670 
km s~^ (Barbon, Cappellaro, & Turatto 1989). 
Adopting the value for Ho of 72 km s^^ Mpc^^ 
consistent with with WMAP and Hubble Space 
Telescope measurements (Spergel et al. 2003; 
Freedman et al. 2001) yields a distance to 
SN 1988Z of 93 Mpc. Previous papers (Fabian & 
Terlevich (1996), Aretxaga et al. (1999), Table 1) 
assumed a smaller Hubble Constant value (Hq = 
50 km s^^ Mpc^^) and thus inferred a larger dis- 
tance (and thus a higher luminosity). 

Table 3 lists the inferred 0.2-2.0 and 0.5-2.0 keV 
fluxes and luminosities for the ROSAT and Chan- 
dra observations. The inferred flux for the ROSAT 
observations is essentially unchanged, as the re- 
duction due to use of the lower Galactic Nh is 
offset by the increase due to the higher kT. The 
luminosities for the ROSAT observations are re- 
duced from those in Aretxaga et al. (1999) by 
about a factor of two as a consequence of the 
revised distance. From the original estimate of 
2±0.7xl0'*i ergs s^^ (Fabian & Terlevich 1996), 
and the adjusted value of 8+3.2x10^"^ ergs s~^ 
(Aretxaga et al. 1999), the peak observed 0.2- 
2.0 keV luminosity is now 3.4±1.2xl0^" ergs s~^. 
Of course, estimates of luminosity in other bands 
also drop by a factor of two as a result of the re- 
vised distance. Despite this downward revision, 
SN 1988Z remains the most luminous supernova 
with a confirmed X-ray detection. 



The 0.2-2.0 keV fluxes are displayed graphically 
in Figure 3, the X-ray light curve. The ACIS flux 
is an order of magnitude lower than the initial 
ROSAT BRl flux. A fit to dlog(flux)/dlog(t) for 
the 0.2-2.0 keV band yields a slope of -2.6±0.6. 
Thus the luminosity evolution of SN 1998Z in this 
band follows the curve Lxoc t-^-^io.e SN1970G, 
with a slope of 2.7±0.9 is the only X-ray super- 
nova with a light curve as steep, but much later in 
its evolution (Immler & Kuntz 2005). 



3.3. Inferred X-ray Properties 

It is possible to infer information about the density 
of the gas from the X-ray luminosity, if an assump- 
tion is made about the size of the X-ray emitting 
region. Aretxaga et al. (1999) show that the Ha 
width decrease is consistent with the shock veloc- 
ity being proportional to t^'"' ". If we assume that 
the Ha width is a measure of the shock velocity, 
it is possible to estimate the radius and hence the 
volume of the region enclosed by the shock for each 
the four X-ray observations. The inferred velocity, 
radius and volume are listed in Table 4. 

For a thermal plasma, Lj, = A(T)ngV2,, where 
P(T) is the plasma emissivity, Ue the electron den- 
sity, and Yx the volume of the emitting material. 
The product n^ Va; is referred to as the volume 
emission measure. For a plasma with a tempera- 
ture in the range 0.5-10.0 keV, the emissivity in 
the 0.2-2.0 keV band is -lO^^ erg cm^ s'^ (Ray- 
mond, Cox, & Smith 1976). In Table 4 we have 
tabulated the value of n^ V^; for each of the obser- 
vations. The emitting volume V^ is equal to V0, 
where V is the total volume enclosed by the shock 
and (j) is the volume fllling factor. Table 4 indi- 
cates that the value of Ugi/)^'^, the density of the 
shock-compressed gas, varies from 6.2x10^ cm^^ 
for the first i?OS'ylT observation to 1.1x10^ cm~^ 
for the Chandra observation. Thus the density 
is high, even for the lowest density case, assum- 
ing the plasma uniformly fills the shocked volume. 
Prior to being shock heated, the density of the am- 
bient gas will be at least a factor of 4 lower than 
the post shock density. 

In Table 4 we also list the inferred X-ray emit- 
ting mass. In all cases, the minimum mass is 
several solar masses. The high shock velocity 
(>1,000 km s~^) ensures that the bulk of the lu- 
minosity emerges in the X-ray band throughout 



the interval covered by the RO SAT and Chandra 
observations. 

The methodology used above for inferring proper- 
ties of the X-ray emitting gas are similar to those 
first used by Imniler and co-workers (as summa- 
rized in Immler & Lewin (2003)). The funda- 
mental difference is that these authors assume a 
constant shock velocity, whereas we assume one 
decreasing with time. If we were to assume a con- 
stant velocity, the values we infer for density and 
mass would be similar, but would be higher or 
lower, depending upon the velocity we assume. 

The rapid decrease of the X-ray emitting mass can 
in principle be used to set a lower limit on the 
density and thus constrain the filling factor. A de- 
cline in the emitting mass indicates that the shock- 
heated material is cooling quickly. The cooling 
time tc can be written as: tc = 0.2 v|/(n7A23) yr, 
where vg is the shock velocity in units of 10^ cm 
s^^, n7 is the density in units of 10^ cni^"^, and 
A23 is the plasma emissivity in units of 10"^"^ 
erg cm^ s~^ (Aretxaga et al. (1999)). The 
contents of Table 4 indicate n7>0.01 and V8>1, 
adopting A23 ^1. That the X-ray emitting mass 
is decreasing at a rate between ~0.2 and ^1 so- 
lar mass per year suggests a cooling time on the 
order of 1-5 years. Such a cooling time requires a 
density >4xl0^ cm^. As this density lower limit 
is consistent with the inferred X-ray emitting gas 
density, we conclude that a filling factor as large 
as unity is allowed. 

It is also possible to estimate the mass loss rate of 
the progenitor, if we assume that the circumstel- 
lar envelope is due to mass loss. Integrating the 
X-ray emitting mass yields a variable mass loss 
rate. For the epoch of the initial ROSAT obser- 
vation, corresponding to ^4000x(vto/10 km s~^) 
years before the explosion, the minumum mass 
loss rate (from assuming (t>=l) is ^1.2x10""^ Mq 
yr~^Vu,/10 km s~^). For the epoch of the Chandra 
observation, which corresponds to ~5500x(v^/10 
km s""'^) years before the explosion, the minu- 
mum mass loss rate is ~4xl0~^ Mq yr~^(v^/10 
km s^^). Alternatively, the X-ray cooling rate 
requires at least 7 M0 to have been swept up. 
This corresponds to an average mass loss rate of 
~1. 3x10-3 M0 yr"i(v^/10 km s^^). This range 
of mass loss rates lie at the extreme end of those 
inferred for X-ray supernovae (Immler & Lewin 
2003), as would be expected given its high lumi- 



nosity. 

We assumed above that the gas is in colli- 
sional ionization equilibrium. For a shock- 
heated plasma, CIE is attained when the prod- 
uct net>10^3cm~3 s. The inferred gas densities 
require t--l-8xlO®0^/^ s = 3-24(/)i/2 yr to reach 
equilibrium. Thus the equilibration time is com- 
parable with the cooling time, and the assumption 
of CIE is justified. 

The X-ray fluxes require Ugi/)^'^ to decline with 
time. Since Lj;Cx t^^'® and Ugi/)^' ^oc L^; , we ex- 
pect Ue^^/^cx t^^-3. The data do not allow de- 
termination of whether the density declines while 
the filling factor remains constant, the filling fac- 
tor declines while the density remains constant, or 
both change with time. In the discussion below, 
we consider the implications of a decline of density 
and filling factor in the context of models proposed 
for SN 1988Z. 



Discussion 



Two distinct models have been proposed to ex- 
plain the multiband light curve of SN 1988Z. One 
postulates that the supernova exploded in an ex- 
tremely dense ('^lO'' cni^^) medium (Fabian & 
Terlevich (1996); Aretxaga et al. (1999)). As 
the supernova shock propagates through such a 
medium, it loses energy rapidly, and evolves in a 
matter of years to the "shell forming," or radiative 
stage. Extremely high luminosities are produced 
as the SN radiates away most of its kinetic energy 
over decades. The second model postulates propa- 
gation of the SN shock through a clumpy medium, 
formed by an unstable wind from a progenitor 
star that underwent extreme mass loss (Chugai & 
Danziger (1994); Wilhams et al. (2002)). In this 
case the forward shock propagates rapidly out- 
ward, driving slower shocks into the clumps it en- 
counters. X-ray (and radio) emission arises from 
the shocked clumps. Below we address how the 
predictions of these two models compare with the 
Chandra observation and the inferred X-ray light 
curve. 

Aretxaga et al. (1999) provide a capsule sum- 
mary of the so-called " compact supernova" model, 
and compare optical, radio, and X-ray data from 
SN 19988Z with the model predictions. They show 



that the hght curves in many bands (for the first 
decade after the explosion) follow trends similar 
to the model predictions. They also show that the 
shock velocity, which they infer from the width of 
the Ha line, follows the t'^/^ behavior predicted 
from the model. (We have used this inference 
above in our estimate of the radius enclosed by the 
shock.) They predict that the bolometric luminos- 
ity should decline as t^^^/^. As their model also 
predicts that the bulk of the luminosity is emitted 
in the X-ray band, the X-ray luminosity should 
follow this decay curve. The ROSAT fluxes ap- 
pear to fit this pattern. They estimate an average 
density on the order of 10^ cm~'^. 

This model runs into difficulty when it is com- 
pared with the extended X-ray light curve. The 
measured t-2-6±o.6 jg ^^^ consistent with the pre- 
diction. The assumption of a constant density of 
~10^ cm^"^ is challenged by the decline in nt/)^/^ 
required by the X-ray flux measurements. A con- 
stant density of ^10 cm^'^ requires a small filling 
factor that decreases with time (from 3.6x10^"^ 
at the time of the first ROSAT observation to 
1.2x10^^ at the time of the Chandra observation). 
If the X-ray emitting gas is confined to a shell, 
its thickness also decreases, from ^2x10^'"' cm to 
~lxlO^^ cm. Additionally, the model is incon- 
sistent with the late-time radio spectrum, which 
shows little free-free absorption (Williams et al. 
2002) from the intervening material composing the 
cocoon in which the supernova exploded. 

There are two ways to potentially reconcile these 
difficulties. First, it is conceivable that the bulk of 
the luminosity has shifted out of the X-ray band. 
This is not likely, based on the inferred shock ve- 
locity from the Ha line width measurements. Even 
if the velocity is overestimated by a factor of two 
at the time of the Chandra measurement, the bulk 
of the emission should still arise in the X-ray band. 
Second, it is possible that the supernova shock 
has encountered a change in the average medium 
density. This possibility cannot be ruled out, 
given the small number of measurements post- 
dating Aretxaga et al. (1999). For instance, it 
is possible that the shock recently broke through 
the dense cocoon into a low-density medium. In 
that case, the X-ray light curve, which until the 
Chandra observation could be characterized by a 
-11/7 slope, might show a pronounced break. If 
the change occurred at day 4500, roughly mid- 
way between the final ROSAT and the Chandra 



observations, then the light curve from day 4500 
to the Chandra observation would have and ex- 
tremely steep decline: Lxoc t~^. The two steepest 
light curves observed are those of this remnant and 
SN 1970G (Immler & Kuntz 2005). Therefore, 
while the Chandra flux measurement does not ex- 
clude the compact supernova model, the required 
steep late-time decline entails significant modifica- 
tion of the model. 

Williams et al. (2002) champion the alterna- 
tive model, whereby the high radio luminosity 
is produced by slow shocks driven into dense 
clumps of wind material. They assume a constant 
shock velocity through the interclump medium 
(20,000 km s^^, as was observed early in the evo- 
lution of SN 1988Z -Turatto et al. (1993)). Using 
this assumption, and the theoretical framework of 
Chugai & Danziger (1994), they interpret the ra- 
dio light curve as arising from a medium with a 
low cloud number density (filling factor). They 
quantify neither the filling factor, nor the density 
of the clumps or interclump medium, but they in- 
fer a mass loss rate of 1.2x10"^ Mo yr~^, and a 
corresponding M /wy^ind = 7.3x10^'^ g cm^^. 

The key observation driving this interpretation is 
that the radio light curve between days 1750 and 
5000 has a -2.7 slope. This slope is essentially 
identical to the X-ray decline between days 2355 
and 5678. While SN radio and X-ray light curves 
are generally correlated, and all X-ray supernovae 
are detected in the radio, none has been observed 
to have exactly the same slope (Immler & Kuntz 
(2005)). In fact, the X-ray light curve should have 
a shallower slope. Since the X-ray and Ha fluxes 
depend on n^ while the radio light curve depends 
on Ug-^, the X-ray light curve should follow the Ha 
light curve, which Williams et al. (2002) show is 
consistent with a -2.1 slope. For our case, since 
the filling factor decreases while the volume in- 
creases with time, a simple interpretation might 
be a shock running through a clumpy medium. 
Based on the clumpy medium model, a light curve 
with slope -2.1 is marginally consistent with the 
X-ray data. 

As a matter of interest, Williams et al. (2002) 
show that the radio and Ha light curves steep- 
ened significantly after day 1750. They interpret 
this steepening as evidence for a decrease in the 
number density of clumps. As none of the X-ray 
observations were performed before day 1750, we 



do not know whether the X-ray Ught curve also 
changed slope. Whatever the physical reason be- 
hind it might be, the fact that one slope change 
has been observed, however, supports the possi- 
bility of a second steepening during the last few 
years. All the observations listed in Williams et 
al. (2002) were performed before day 4500. If 
we assume that the X-ray light curve followed a 
-2.1 slope until that time (allowed by the ROSAT 
fluxes) , and then decreased more rapidly after that 
date, then a slope of at least -2.8 is required for 
consistency with the Chandra flux point, which 
corresponds to a radio slope of -3.9. Again, no X- 
ray light curve this steep has ever been measured. 

In summary, neither model is consistent with all of 
the observations. Aretxaga et al. (1999) assume 
a dense, uniform shell. This should result in self- 
absorption of the radio emission which becomes 
smaller with time. Williams et al. (2002) show 
that the self-absorption is small, and constant, 
suggesting radio emission is not being produced 
within a dense cocoon. If the intense UV and X- 
ray radiation keep this material ionized, then little 
or no intrinsic absorption of soft X-rays might be 
detected at this time. The highest allowed col- 
umn density from the quantile analysis is 5x10^^ 
of neutral material. Williams et al. (2002), on the 
other hand, assume a constant, high shock veloc- 
ity of 20,000 km s^^. This is not consistent with 
the smooth velocity decline suggested by the Ha 
width. From the X-ray standpoint, the light curve 
is at odds with both models. Both predict flatter 
X-ray decays than observed: -11/7, and -2.1 com- 
pared with the measured -2.6. But the current 
X-ray observations do not rule out either model. 

The ambiguity presented by the limited X-ray 
data could be largely removed by a single addi- 
tional observation. Such an observation will show 
whether the light curve has a shallower average 
slope than inferred from the existing X-ray obser- 
vations, continues to track the -2.6 slope, or steep- 
ened after the i? 05"^ T" observations. If the slope 
is flatter, then it could support either model and 
possibly provide a discriminator between them, 
depending on whether it is closer to the -11/7 
predicted by the compact supernova model or - 
2.1 predicted by the clump model. As described 
above, a change in slope could rescue either model, 
but would require a change in the properties of the 
medium being encountered by the shock. If the 
slope is constant at -2.6, then a different model 



must be sought which explains the X-ray light 
curve in the context of the numerous other ob- 
servations. Additionally, radio observations are 
needed to show whether the X-ray and radio light 
curves continue to follow each other. 

It is possible to compare the inferred density pro- 
file with those of other X-ray emitting SNe. Imm- 
ler et al. (2005) plot the density profiles for sev- 
eral X-ray SNe, and show'^ that ncsM oc r~^, 
where s ranges between 1 and 2. The inferred 
density (or more precisely nxcj)^^"^) follows a r~^'^ 
profile, entirely within this range. In calculating 
their density profiles, Immler et al. (2005) as- 
sume a constant shock velocity. In contrast to this 
assumption, we have assumed that the shock ve- 
locity profile of SN 1988Z is known, and that the 
shock decelerates. A comparison with the inferred 
density for, e.g., SN 1979C indicates that its in- 
ferred density range is comparable to the range of 
nx0i/2 foj. gj^ 1988Z. This could be interpreted 
as suggesting that the actual density is substan- 
tially higher, and thus (f> is much less than unity 
(or that the shock velocity decreases in the other 
SNe, leading to different density profiles). 



Conclusions 



An ACIS observation yields a clear detection of 
SN 1988Z 16 years after discovery, and a decade 
after the first X-ray detection. Although limited 
spectral information is available, the X-ray light 
curve constructed using this observation and the 
three archival ROSAT observations yields a flux 
decline proportional to t"^-^^"-^. This decline is 
steeper than predicted by either of the two promi- 
nent competing models. The inferred density is 
high, larger than 10^ cm"^, indicating a mas- 
sive envelope around the progenitor, and requir- 
ing rapid evolution of the X-ray emitting plasma. 
If the circumstellar envelope is the result of pre- 
supernova mass loss, a very high mass loss rate is 
required, as might be expected from this, the most 
X-ray luminous supernova. An additional observa- 
tion can potentially discriminate between the two 
competing models, or require a new model. 
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•^Note that in the cited paper, the density is represented by 
p; we made the change to conform to our notation. 
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Table 1: Summary of Distance, Column Density, and Model Assumptions 







D 


Nh 




Inferred Lx** 


Authors 


Ho 


Mpc 


cm"^ 


Model 


erg s"' 


Fabian & Terlcvich (1996) 


50 


131 


3x10^'' 


1 keV bremss 
5 keV bremss 


IxlO"*^ 
2x10^1 


Aretxaga et al. (1999) 


50 


131 


2.5x10^0 


1 keV bremss 


SxlO"'" 


This paper 


72 


93 


1.4x1020 


5 keV MeKaL 


3.4x10''" 



^For the bandpass of 0.2-2 keV. 



Table 3: 0.2-2 and 0.5-2 keV Fluxes for SN 1988Z 



Detector Obs. Date^ 



Count 
MJD Age'' ExpT Rate'^ 



Flux'' Luminosity'' 

0.2-2.0 0.5-2.0 0.2-2.0 



HRI 1995 May 16-25^ 49884 2335 12287 

HRI 1996 Dec 14^ 50725 2924 6739 

HRI 1997 May 13-24^ 50908 3085 34300 

ACIS 2004 Jun 29 50908 5678 44618 



11±4 3.92 2.75 3.4 

5±4 1.78 1.25 1.5 

3.7±2 1.32 0.92 1.15 

6.7±1.5 0.36 0.25 0.27 



"MJD at center of observation when spanning multiple days; those observations are footnoted, 
''units = days; based upon adopted date of maximum = 1988 Dec 12 = MJD 47508. 



^Units are 10""' scc^l for 0.2-2.4 keV band. 



Adopted model is an absorbed vmekal, kT 



5.0 



keV, Njj = 1.4x10-^" cm ^ . For the HRI values, fluxes were computed from the model after adjusting the model normalization 
to match the observed count rate, 
■luminosity units are lO**" erg s"'^ for 0.2-2.0 keV band. 



Table 4: Inferred Emission Properties of SN 1988Z 



(days) 


Velocity'' 
(km s-') 


Radius'' 
(pc) 


n^^V'' 
(10^3 cm-3) 


n0i/2 

(10'^ cm-3) 


07^ M, 


2335 
2924 
3085 

5678 


1900 
1620 
1560 
1010 


0.044 
0.046 
0.047 
0.056 


3.4 
1.4 
1.15 
0.27 


5.7 
3.4 
3.0 
1.1 


3.2x10-3 4.9 
1.2x10-3 3.7 
9.1x10-'' 3.2 
1.2x10-4 2.1 















"MJD at center of observation when spanning multiple days; those observations are footnoted. 
''Estimated using velocity profile in Aretxaga et al. (1999). 
''(f) is the dimonsionless volume filling factor. 

'^Volume filling factor assuming n=10'^ cm"^ (Aretxaga et al. 1999). The volume filling factor changes by a factor of 100 for a 
change in the assumed number density by a factor of 10. 



6. Figures 



Fig. 1. — X-ray positions of detected sources over- 
laid on an blue optical Digital Sky Survey 2 image 
of the field surrounding SN 1988Z. SN 1988Z is 
at the center of the image; its host galaxy, MCG 
-1-03-28-022, which is not an X-ray source, lies due 
west. 
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Fig. 2.— Quantile plot for SN 1988Z using the 
description of Hong, Schlegel, & Grindlay (2004). 
The underlying grid represents a MeKaL model of 
a coronal equilibrium plasma; the 'mountainous' 
shape arises from changing line emission as (pri- 
marily) the temperature parameter is varied. The 
Nh axis starts at 0.001 x 10^^ cm~^ with isolines at 
0.01, 0.1, 0.5, 1.0, and S.OxlO^^ cm'^. The tem- 
perature axis shows isothermal lines of 0.2, 0.4, 1, 
2, 5, and 10 keV. The position of SN 1988Z cor- 
responds approximately to Nh ~0.5x10^^ cm~^ 
and kT ~5 keV. The temperature range includes 
~3-10 keV; the Nr range is ~0.01-2xl022 cm'^. 



Table 2: ACIS Count Rates for SN 1988Z in vari- 



ous bands 

Band 


Net 
Counts 


Rate^ 


0.2-2.0 


29.8±6.82 


6.68±1.53 


0.3-2.0 


30.4±6.82 


6.81±1.53 


0.3-10.0 


37.4±7.71 


8.38±1.73 


0.5.2.0 


29.6±6.60 


6.64±1.49 


2.0-10.0 


7.7±5.13 


1.73±1.15 


Units = 10~* counts 


s-i. 





o 




Nh-0. 001 .0.01 



log m 



10 



Fig. 3. — X-ray light curve in the 0.2-2 keV band 
based upon ROSAT RRl and Chandra ACIS ob- 
servations. The data points have all been con- 
verted to flux using the same model. The best 
fit power law yields a light curve fcxt~^'^. Also 
shown is a power law slope of -11/7, as predicted 
by Aretxaga et al. (1999). 
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